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Abstract 

The possible boundary conditions consistent with the integrability of the classical sine-Gordon equation 
are studied. A boundary value problem on the half-line a; < with local boundary condition at the origin 
is considered. The most general form of this boundary condition is found such that the problem be 
integrable. For the resulting system an infinite number of involutive integrals of motion exist. These 
integrals are calculated and one is identified as the Hamiltonian. The results found agree with some recent 



1/^ ' work of Ghoshal and Zamolodchikov. 



> ■ 1 Introduction 

The problem of factorizable scattering on a half-line necessarily involves the so-called "reflection equation" 
H, ^ . This equation supplements the Yang-Baxter relation for the factorization of three particle scattering 
and represents the factorization of two particle scattering in the presence of a one-dimensional 'wall' at a; = 0. 
I Recently there has been increased interest in this problem, and solutions have been suggested for sine-Gordon, 
a^ '• afiine Toda, and 0{N) sigma models |, |, |], ||. 

Equations of "reflection" type are also used within the framework of the quantum inverse scattering method 
to prove the integrability of open quantum spin chains, i.e. chains without periodic boundary conditions 
Q^l I ^ I ■ The classical limit of these results was developed in a paper by Sklyanin |jri| , where it was found 
1^ ■ that they generalize the Hamiltonian, (r-matrix), approach to inverse scattering theory jl^, to the case of 
nonlinear partial differential equations with local boundary conditions. The adaptation of these results to 
inverse scattering theory on the half line x < with a local boundary condition at x = 14 saw 
the problem continued evenly to the whole line with a 'point-spin' impurity at the orgin. These results 
are equivalent to the approach in |l^, [l^ where symmetrical reductions of Backlund transformations were 



X 

H , ir-\ 

I considered. The analysis of these problems on a finite interval is to be found in [17, 

Ghoshal and Zamolodchikov j^, conjectured both the classical and quantum integrability of the sine- 
Gordon theory on the half line with a boundary term at the origin. A classical action is. 



S — dt dx 

J —oo J — oo 

the minimisation of which leads to the boundary value problem. 



1 2 

0(-a;) ( (l-cos/30(a;,i)) ) + <5(x) t)) 



(1) 



d''d^(f>{x, t) + — sin /30(x, t) = x < 0, 



f3 

d(j){x,t) dV{(f>{x,t)) 



= a; = 0, 



dx d4> 

lim (mod27r//3)) = 0, (2) 



X — > — CXD 



with the boundary value at — oo being attained in the sense of Schwartz, i.e. (j) ~ 0(mod27r//3) and all its 
derivatives decay faster than any inverse power of a; as x — > — oo. 
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By demanding that a modification of the first non-trivial integral of motion of the bulk theory remain 
conserved on including the boundary it was found Q that in general, 

mx,0) = ^cos(M^4^), (3) 

where A and 00 £^re arbitrary real constants. In |^ the existence of an infinite number of integrals of motion 
for the system (||),(H) was assumed and the associated quantum field theory was studied. 

In this paper, the question of integrability of the classical system (|^) is addressed. The results of Sklyanin 
are used to prove that an infinite number of integrals of motion do exist, but only for certain V^cf)). It is found 
that (||) is the most general boundary term compatible with the integrability of (||). The results thus agree 
with Ghoshal and Zamolodchikov. Further analysis of the system @,(||) can be found in JTof . 

In section ^ the auxiliary linear problem of inverse scattering theory is explained. The case of periodic 
boundary conditions is considered first, followed by the calculation of the integrals of motion for sine-Gordon. 
Next follows Sklyanin's extension to other types of boundary value problem. 

Finally, in section ^ the formalism developed is applied to the system (^) , and the requirement on the form 
of the boundary term V{4>) is deduced. The integrals of the resulting boundary value problem are studied and 
the Hamiltonian is identified. 

The notation of |13] is used throughout the following sections. 



2 The auxiliary linear problem of inverse scattering theory 
2.1 Periodic boundary conditions 

The zero-curvature approach to inverse scattering relies on the existence of a pair of linear partial dif- 
ferential equations for a, d x d matrix whose elements are maps from R^xC to C. The solution of such a 
system subject to a certain initial condition defines the transition matrix T(x,y,t, a) uniquely. Suppressing 
the argument of T{x,y,t,a), and understanding that x >y throughout, the pair of equations and the initial 
condition are. 



dT 

dx 
dT 

'dt 

T{y,y,t,a) = I. (6) 



U{x,t,a)T (4) 
V{x,t,a) T ~T V{y,t,a) (5) 



Here, U{x^t,a),V {x,t^a) are d x d matrices whose elements are functions of a complex valued field ^(x,t), 
its derivatives, and a spectral parameter a G C. 
The formal solution to (^),(||) is given by. 



T{x,y,t, a) — V exp ^^J U{x' ,t,a) dx'^ 



(7) 



where V denotes path ordering of noncommuting factors, jl^ . 
Compatibility of (|) and (||) gives the zero curvature condition, 

dU dV 

----,[C/,y]=0, V.EC (8) 

and this should imply the desired nonlinear evolution equation for $. 

For the sine-Gordon equation, U and V are 2x2 matrices and can be written in terms of the standard 
Pauli matrices crfc, k = 1, 2, 3 and a real valued field as, 

, ^ P d(j) m . ( m ( (3(t)\ 

U(x,t,a) = — ^CT3 + 7r<^o^ii"^i'^ ~?r '^i + TT^i^^^'^o^ Mr (9j 



Ai dt 2i \ 2 J 2i \ 2 

V[x,t,a) = — — 0-3 -1- — smha sm I — I CTi + — cosha cos I — ) 0-2- (10) 
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The imposition of periodic boundary conditions on (j){x, t) with period 2L ensures periodicity of V{x, t, a) 
also, and (||) becomes, 

dT{L,-L,t,a) ^ j^^^^ ^^^^^^^ ^^^^ 

Taking the trace of (|ri|) and using the chain rule for differentiation gives, 

d ln(trT(£,-L,t,a)) 

0^ = 0, Va, (12) 

which can be expanded as a Laurent series in a to give an infinite number of integrals of motion, related to 
the field by (0) . The existence of such an infinite set is crucial for the integrability of the original evolution 
equation, in particular for the construction of action-angle variables for the theory. 

In the limit L — > oo, eq.(ll) simplifies as (j) and its derivatives vanish, and the evolution of the matrix 
elements of T, the transition coefficients, have a simple form. A detailed treatment of the auxiliary linear 
problem (Q), and its inversion, the Riemann problem, can be found in p^ . 

The alternative formulation to the one outlined above is the so-called Hamiltonian, or r matrix, approach. 
The problem is now regarded as an infinite dimensional Hamiltonian system with phase space M. Points in 
this space are determined by pairs of complex 2L spatially periodic functions at fixed time ($(a;, i), n(a;, t)), 
n being the conjugate momentum to The 'algebra of observables' ,4 on is defined over C and comprises 
functionals A : M. — >C, bounded everywhere and having smooth functional derivatives that are 2L periodic 
functions of x. The Poisson bracket on this algebra is defined at constant time slices by. 



/M[$(w,t),n(M,i)] 5B[<^{y,t),I{{v,t)] 5A[<^{u,t),IV{u,t)] 5B['^{v,t),IV{v,t)] 



(5n(z,i) 5<^{z,t) (5$(z,i) 5W{z,t) 



dz, (13) 



where A, B E A and u,v,z now lie inside the fundamental domain [— L, L]. This Poisson bracket can be 
generalised to act on d x d matrices, C, D whose elements are in A, by replacing A,B in (|l|) by C^^\D^^\ 
respectively, where C^^-* := C(iil, Z?^^-* := I®D^ and / is the dxd identity matrix. As a result of this notation, 
{C'^^\ D^^^} is a X matrix whose elements are the Poisson brackets of all possible combinations of matrix 
elements of C, £>. 

The Hamiltonian formulation of the problem relies on the existence of an r matrix. This is a function from 
C to End (C^® C^). It is defined in terms of the transition matrix, T{x,y,t, a) ^) by the equation, 

{T(i)(a;,y,t,a),r(2)(a;,y,t,7)} = [n ^(a - 7), T^) (a:, y, t, a)T(2) (x, y, i, 7)] , Va,7eC. (14) 

The dependence of r on the difference {a — 7) is not the case for all models, but it is so for the sine-Gordon 
equation considered below. Antisymmetry and the Jacobi identity for the bracket hold if and only if, 

ri2(a - 7) -f r2i(7 - a) = 0, (15) 
[ri 2 (a - 7)7^13(0;) -I- r2 3(7)] + ['"13 (a), ''23(7)] = 0) (16) 

where the subscripts denote the nontrivial action of r in C'^ C*. Eq([l^) is the classical Yang-Baxter 
equation. 

For sine-Gordon, phase space Msg is the subset of M given by pairs of real 2L periodic functions 

{(f>{x,t),7T{x,t)), 

2tt dd(x t) 

(f>(x + 2L,t) = (l)(x,t) mod —, Tr( x + 2L,t) ^ n(x,t), n(x,t) :^ ' \ (17) 

p ot 

On this space will be defined the algebras Asg and A^^,. The algebra Asg is the restriction of ^ to A^sg, 
while „4sg is the subset of Asg mapping J^sg which becomes an algebra on restricting the scalars of 

Asg to R. For convenience, the Poisson bracket will be defined on both these algebras but it is understood 
that the integrals of motion should be elements of A^^^. Also, d = 2 and, 

rsG{a - 7) 77- — TT ^{1 ® / - 0-3 (g) (73) - . -(di (g) (71 + (72 ® (72). (18) 

Id tanh(Q! — 7) lb smh(a — 7) 
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Returning to the integrals of motion of a general theory, taking the trace of (^) and using the chain rule 
for functional differentiation gives, 

{ln(trT(a;,?;,i,a)) ,ln(trT(x,y,i,7))} = 0, Va,7eC, x,y s.t. - L < y < x < L (19) 

which, on expansion as Laurent series in a, 7 about the singularities in T, gives an infinite number of involutive 
charges, depending on $ through (Q). However, so far no mention has been made of the boundary conditions 
on These have been used in defining the 'algebra of observables' A, and the fundamental domain 
—L < y < X < L. For these involutive quantities to have smooth functional derivatives and so be elements 
of this algebra requires that x — L,y ^ —L. Hence, the integrals of motion found using the zero-curvature 
formulation are also in involution with respect to the Poisson bracket ( |l3|). This suggests that the model may 
be completely integrable, along the lines of the Liouville- Arnold theorem |2C|| for finite dimensional Hamiltonian 
systems. 

Alternatively, there is no need to appeal to the zero-curvature condition at all. By requiring the involutive 
charges to be elements of A, they are fixed uniquely. Among these should exist a Hamiltonian H{L, —L,t) 
such that, 

^5^ = {ff(L,-L,i),<&(x,t)}, ^^^^ = {H{L,-L,t),Il{^,t)}, ~L<x<L, (20) 

imply the field equation. The problem can then be extended to the x axis by letting L — > cx) and demanding 
that <i> be a Schwartz function. 

For the sine-Gordon equation such a Hamiltonian is found as follows. Fixing Ima = n/2, (^) shows 
T(L, — L, t, a) to have essential singularities at Re a — ±00. As a result, it is expected that there will exist two 
distinct infinities of involutive integrals coming from the coefficients of the Laurent expansions about these 
two "points" . First consider Re a — > +00. Suppressing t and defining T'{L, — L, A) = T{L, — i, a) and A := e", 
T'{x, y, A) is rewritten as, 

T'{x, y. A) = ^l{x)f{x, y, X)n-'{y), (21) 
where il{x) = exp{i(3(j){x)a3 / i) . It can be shown ||l^ that for A ±ioo, 

f{x, y, A) = (/ + W{x, A)) exp Z(x, y, A) (/ + Wiy, X))-' + 0(|Ar°°), (22) 



where, 



m-^A) ^ (23) 



n=0 



Z{x,y,X) - ^3+^2^ ^„ • (24) 



n=l 



The Wn , Zn have the form. 



with, 



wo{x)=i, (26) 
2 dw^jx) P ( I dHxy\ . ^ . « f . f s 

^ ^ k=l 



' ^^,(^)^„_,_^(^)_^e^WW5„,i, (27) 



2 

fc=0 



zn{x,y) = — I {w„+i{x')-e-'^^^^'^w„^,{x'))dx'. (28) 
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Now imposing periodic boundary conditions at ±L simplifies the form of tr T'{L, — L, A). As a result, the 
first set of integrals will be the coefficients of tr exp Z(L, —L, A). The unimodularity of T'{L, ~L, A) implies 
that Z{L, —L, X) be traceless up to 0(|A|~°°) terms, and so the z„(L,— L)'s are real and Z{L,—L,X) is 
proportional to (J3. Thus, 

trr'(L,-L,A) = icos(^£f^fc^-^^+0(|^ A ^ to. (29) 

Rewriting the cosine in terms of exponentials and neglecting the one which is 0(|A|~°°), the result is, 

mXL zi{L,—L) Z2{L,—L) 



ln(trr'(X,-L, A)) i 



X A2 



o(iAr 



A +icx), (30) 



and this series generates the set of charges {z„(L,— L)} defined by (|27|),(P8D. These charges are purely 
imaginary, time independent, in involution with each other, and are elements of Asg- To construct elements 
of it suffices to multiply all the charges by i. Note that (|2^) , (^) also hold for a — — 7r/2 or A — > —ioo. 
The resulting set of charges in Asg are complex conjugate to those found above and lead to the same elements 
of on multiplication by —i. 

The second distinct set of integrals comes from the expansion of tr T'{L^ — L, A) as A ^ iO±. To find this, use 
is made of the invariance of (^),(0) under the replacement 7r(x,t) — > 7r(a;,i), (l){x,t) —<j){x,t), e" — e~". 
As a result 

f ' 



T' i^x,y,-jj =T'{x,y,X), (31) 

where T'{x,y,X) denotes the transition matrix for 7r(a;,i) = Tr(x,t), (j){x,t) = —(f>[x,t). So taking A — s- «0± in 
(31) amounts to taking ^ := — 1/A ±ioo in T'{L, — L,/i). Therefore the analysis given above for A — > ±ioo 
can be used again on the replacement A —1/X,4>{x,t) -~4>{x,t). The resulting conserved charges are, 

z„(L, -L) = ^ J\w^+i{x') - e^^^^^'^w^^,{x'))dx', (32) 

with the WnS given by (|27| ) on replacing (j) —(f). Repeating the same procedure as above for this second set 
shows the {zn{L, —L)} to be elements of -4*^. 

The linear combination H{L, — L) := ^(zi(L, — L) + zi(L, —L)) is of the form. 



(1 -cos /?(/)) dz, (33) 



and is identified as the Hamiltonian of the model. The sine-Gordon equation results from (^0|). Note that all 
the Zn{L, — L), Zn{L, —L) are zero for n even. Also note that if is an even function of x then the Hamiltonian 
density is itself even. In general the densities Zn{L, —L) + Zn{L, —L), n = 1, 3, .. have this property, ("parity 
even"). Finally extending L —> 00 and demanding that (j) have Schwartz boundary conditions, gives the theory 
on the whole line. 



2.2 Local boundary conditions 

The analysis given below was first provided by Sklyanin . The idea is to modify the Hamiltonian approach 
of the last subsection by introducing a new 'algebra of observables' B, the form of which is not specified. This 
is to be contrasted with the precise definition of the algebra A given previously. By analogy with the transition 
matrix T(x,y,t,a), a, d x d matrix with elements in B is introduced and called T(x,y,t,a). This matrix is 
required to satisfy the equation, cf . (p^ , 



{r(i)(a),r(2)(7)} = [ri2(a-7),rW(a)r(2)(^) 



T(i) (a)ri 2(a + j)T^^^ (7) - T^^) {j)r, 2 (a + ^)T^^'> (a), 

(34) 
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Vq;,7, and the x,y,t dependence has been suppressed. The Poisson bracket is defined in (|T^) with the 
exception that the range of the integral is changed to [x_,a;_|_], and it is understood that < y < x < x-^-. 
Antisymmetry of this bracket implies, 

ri2{a) - r2i{a) = 
ri2(-a)+ri2(a) = 0. (35) 
Constant representations of the Poisson bracket algebra B must also satisfy ( p^ , hence 

ri2(a-7),^i'^(«)i^f (7)] +/4'H«K2(« + 7)^i'^(7)-^f (7)^i2(« + 7)i^i'^(«) -0, (36) 

where K± (a) are d x d C numerical matrices of the complex spectral parameter a and possibly some complex 
constants. From such a K-[a), a T[x,y,t,a) satisfying (^) can be built according to, 

T{x,y,t,a) = T{x,y,t,a)K_{a)T-^{x,y,t,~a), (37) 

where T{x, y, t, a) is given by (0). 

Multipying T{x,y,t,a) on the left by K+{a), evaluating the bracket of this with itself using (|3^), taking 
the trace and using the chain rule for functional differentiation gives, 

{hi{tr K+{a)T{x,y,t, a)) ,hi{tr K^{j)T{x,y,t,j))} ~ 0, Va,7GC, x,y s.t X-^<y<x<x+, (38) 

i.e. an infinite set of involutive quantities. The question remains, however, as to whether these can be used to 
construct elements of the unspecified algebra B. As was the case for the algebra A, the requirement of smooth 
functional derivatives of elements of B imposes y = x = In principle, specification of the boundary 
conditions at x± on elements of phase space should be sufficient to determine B uniquely. In contrast to the 
periodic case, however, local boundary conditions of the form + dg>V{^)\^^^^ = do not impose any 
general form for functional derivatives of elements of B. So, given boundary conditions at x±, the entire form 
of B cannot be deduced on symmetry grounds, there are no constraints on the forms of the K±(a), and a 
construction analogous to the one leading from ( p^ ) to (^) cannot be performed. 

The way around this problem is to appeal to the zero-curvature formulation. For periodic boundary 
conditions, ( p^ can be regarded as a consistency check on the form of tTT{x,y,t, a) such that it be an 
element of A. In this case, such a check just imposes y — ~L,x — L which gives no new information on the 
structure of the algebra. If the same point of view is taken for the nonperiodic case, however, K±(a) become 
constrained by possible forms of boundary condition at x±. To see this, consider, dtK^(a)T{x^, x-,t, a) 
using (^),(^7|). The trace of this quantity should be zero if the Hamiltonian and zero-curvature approaches 
are to be consistent. This requirement leads to, 

K±{a)V{x±,t,±a) ^V{x±,t,Ta)K±{a), Va G C, (39) 

which contains nontrivial information as to the forms of the K±(a), for particular boundary conditions of the 
theory. 



3 Application to sine-Gordon theory 



In this section the formalism developed in subsection |2.2| is applied to the sine-Gordon equation. For sine- 
Gordon d = 2 and the matrices U,V, appearing in the zero-curvature condition are given by (p|),(p^. The 
basis of EndC^ is chosen to be the Pauli matrices and the two dimensional identity ctq), and K±(a) is 
written in this basis as Z?^^(a)cr^, = 0, ..,4 . It should be recalled that the standard Pauli matrices satisfy 



aiCTj = ieijkCTk + 5ij<J{), fc = 1, .., 3. Substituting K±, V into ( p9D gives. 



/3 



2i 



/30 



^±(") ( -^-Q^ [f^M'^s] + — sinhasin ( — ) {cr^, fii} -I- — coshacos ( — ) [cr^,cr2] 



2i 



0. 



The anticommutator implies D}j_{a) = 0, and so the most general form this equation can take is, 



9^ 
dx 



2mP±{a) 



PI 
2 



2mQ±(a) 



PI 
2 



(40) 



(41) 
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where, 



Diia) 



Q±{a) 



Dlia) 



cosh a. 



Alternatively, K±{a) = D'^ao implies. 



with. 



i?±(a)m 
2/3 



sm 



2 



= 0, 



R±{a) 



Dj_{a) sinh a 



(42) 



(43) 



(44) 



From now on only the boundary at a;-)_ will be considered; however similar arguments would apply to the 
other boundary. The subscripts ± can therefore be suppressed. Eq.(pl|) would be in the form of a boundary 
condition that could be imposed on (p if P{a), Q{a), were real and independent of a. Similarly (|4^) would be 
a boundary condition if R{a) were independent of a. Denoting these constants by P,Q,R respectively, the 
form of K{a) is constrained even more. The result is. 



K{a) = - 



D\a) 



sinh ( 



(Pctq + i sinh aa2 + iQ tanh atra) , 



or, 



K{a) 



iR 
sinh c 



-(Jo- 



lt can be verified that these satisfy ( p^ ) so that ( |38| ) holds with (p5[) , (|4^) . 

Since the K{a) are only defined up to an overall scale, it will be convenient to work with, 



k{a) 



P 



sinh a 



(To + 1(^2 



cosh a 



0-3, 



(45) 
(46) 

(47) 



and 



The normalization is fixed by requiring det K{a 
boundary conditions of the form. 



K{a):^ao. (48) 
1 as Re a — > ±cxd. So, by demanding that (j) satisfy 



2m dx 



or, 



P± sin 
2n7r 



/30 



Q± cos 



0, 



-0, 



n = 1,2,, 



(49) 



(50) 



the associated K±{a) can be found and an infinite number of integrals in involution can be constructed. These 
will be the coefficients in the expansion of, 

ln(^trX+(ln(A))T'(a;+,a;_,A)X_(ln(A))T'"\a:+,2;_,l/A)) , (51) 

as A ^ ibicxo, iO±. From these charges, elements of the algebra of observables S*q will be constructed as they 
were for periodic boundary conditions. 

Now, consider the boundary value problem (H). Taking — — oo, (|5^) becomes the required Schwartz 
boundary condition. With a;+ = 0, ( ^9| ) has the form of (||),(^) after relabelling constants, and the desired 
boundary value problem has been fitted into the developed formalism. Hence the conclusion that (||), (||) is 
integrable. It remains to construct the algebra of observables and identify a Hamiltonian for the problem from 
the infinity of integrals. This is done below. 
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First put the problem (^,(|^) on the finite interval [— L,0] and define 

P iQ 

Ka{a) := — — (Tq + i(^2 H -, — era, (52) 

smh a cosh a 

K^L{a):=ao, (53) 

so that dsi] ) becomes, 

In (tri^o(ln(A))T'(0, -L, \)k^LT'{-L, 0, 1/A)) . (54) 
To expand this in an asymptotic series as A ^ ±ioo, use is made of (p]) , (|2l|) , (p2[) giving, 

In (tri'o(ln(A))r2(0)(/ + W{Q, A)) exp Z(0, -L, A)(/ + W{-L, X))-^n-\-L)n{-L) 

x{I + W{-L,-X))expZ{~L,0,-X){I + W{0,~X))-'^n-\0) + O{\X\-"°)^ , (55) 



where, W, Z are defined by (|23|)-(^8|) and the hats indicate the same quantities but with (f> ~^ —(/). Imposing 
the Schwartz boundary condition at — L simplifies this further, the result being. 

In (tri^o(ln(A))T'(0, -L, A)/C_L(ln(A))r'"'(0, -L, 1/A)) 

= In (tr (/ + W{0, -A))-il7(0)i^o(ln(A))f7(0)(/ + W{0, A)) exp(Z(0, -L, A) - Z{0, -L, -A))) (56) 

mod(0|A|"°°) as A ^ ±ioo. From (||), it follows that, 

Z(0, -L, A) - m L, -X)^^± ^^^^^a, - ± ^^^^7, (57) 

n— — 1 n— 1 

where 2;„(0, -L) - (-l)"z„(0, -L) = a„(0, -L) + i5„(0, -L), n = 1, 2, .., and ao(0, -i) = 0, 

a-i(0, -L) = -mAL/2. Hence, only the diagonal terms in (/ + W{0, -X)y^n{0)Ko{\n(X))n{0){I + W{0, A)) 

will contribute to the trace. Defining, 

diag ((/ + W{0, -X))-'mKo{HX)MO)iI + W{0, A)) = £ (^^3 + ^l) , (58) 



where (p2),(p6[) and the involutions. 



(/ + w{o, -x))-^n{o)Ko{HX))n{o){i + w{o, x)) 

= (T2{I + W{0, -Xj)-^n{0)ko{\n{X))n{0){I + W{0, X))a2, (59) 



exp(Z(0, -L, A)) - Z{0, -L, -A) = (T2 cxp(Z(0, -L, A) - Z(0, -L, -Xj)a2, (60) 
are used to show that the {c„(0), rf,i(0)} are real. Thus (^) has the form, 

\ n=0 \,n=-l I 11=0 \m=-l / / k=l 

mod(0(|A|^°°) as A ^ ±ioo, the Laurent expansion of which has real coefficients. Now consider the expansion 
of this as A ^ +ioo. Rewriting the sin, cos in terms of exponentials and discarding the ones which are 
0(|A|-°°), ^) becomes, 

* 2^ Xm 2^ ;^/c + 2^ I ' 

m=-l fe=l \n=0 / 



as A ^ +ioo. From (|58|),(52),(|23|),(P5[) one finds cq = 1, do = 0, so that the logarithm can be expanded as 



a power series. The coefficient of (p3) at 0(1/A) is i{ai ~ di) — (61 — ci). This is, in principle, a complex 
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valued functional and is defined to be an element of Bsc- Indeed, this will be the case to all orders. It remains 
to construct the integrals of motion which are elements of B^^. To do this, use is made of the A — > —ioo 
expansion of (^) which has coefficients at each order in A that are complex conjugate to those in the A —> +ioo 
expansion. Hence, not only are the complex integrals in involution with one another, but also the complex 
conjugated ones too. (This is an obvious consequence of the real coordinates on M-sg)- As a result, the real 
and imaginary part of each integral is separately in involution with the real and imaginary parts of all the 
other integrals. Therefore once a Hamiltonian for the system is identified as being the real or imaginary part 
at some order in A, then the real or imaginary parts of the others can be chosen at will. Thus, by multiplying 
terms by i an infinite number of elements of in involution, can be constructed. At 0(1/A) in (|6^), the 
first element of Bsc is, 



-1(3^ 
2m 



1 n 1 fd(l)V , 



/30(O) 4Qm . /30(O)\ 
— sm — - — 



4Frn f-^w-y r^-ry^, 1 / N 

'^^+^cos^ -p^^^^^-j- (63) 



Multiplication of this by 2mi/ 0^ defines an element of and is identified as the Hamiltonian of the system 
(||),(^ defined over the interval [— L, 0]. By imposing Schwartz boundary conditions at —L the integral can 
be extended to the whole half-axis. 

The 0(1/A^) term is found to be real and equal to 4(P^ + Q^), i.e. a trivial integral of motion. The 
coefficient at 0(1/A'^) is purely imaginary and is made up of the integral over [—L, 0] of the next "parity even" 
density plus a boundary term. It is not clear as to whether the series continues to alternate between real and 



imaginary terms. However, it will be the case that the integrals at 0(1/ A 



2n\ 



1,2, .., will depend purely 



on the boundary at a; = 0. 

Finally it remains to return to ( |5^ ) to calculate its expansion as A ^ «0±. Replacing A by — in ( ^ ) 
and using ( |3l| ) one finds that the analysis for A ±zoo will hold on replacing (j) — > — ^, Q ~Q. So all the 
integrals previously found will form another infinite set on imposition of this transformation. However the 
first three are invariant under this, and hence it is likely to be the case in general. 



4 Conclusions 

The system (§),(||) conjectured to be integrable by Ghoshal and Zamolodchikov has been fitted in to a formal- 
ism developed by Sklyanin for integrable equations with local boundary conditions. It has been shown that 
an infinite number of integrals of motion, in involution, do exist. These charges fall into two classes: those 
consisting of integrals of "even parity" densities plus a boundary term, and those depending purely on the 
field and its derivatives at the boundary. 

This analysis fits into the scheme developed by KhabibuUin on symmetrical reductions of Backlund trans- 
formations, but it remains to continue the problem evenly to the whole axis with a "point-spin" at the origin 
along the lines of . The analysis of the solutions of such a system is an interesting open problem. 
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